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Abstract 

The problem of stability of saturated and non-saturated ferromagnetism in the 
Hubbard model is considered in terms of the one-particle Green's functions. Ap- 
proximations by Edwards and Hertz and some versions of the self-consistent approx- 
imations based on the 1/z-expansion are considered. The account of longitudinal 
fluctuations turns out to be essential for description of the non-saturated state. The 
corresponding pictures of density of states are obtained. "Kondo" density-of-states 
singularities owing to spin-flip processes are analyzed. The critical electron concen- 
trations for instabilities of saturated ferromagnetism and paramagnetic state are 
calculated for various lattices. Drawbacks of various approximations are discussed. 
A comparison with the results of previous works is performed. 

71.10.Fd Lattice fermion models (Hubbard model, etc.), 71.28.+d Narrow-band 
systems; inter mediate- valence solids, 

1 Introduction 

The problem of ferromagnetic ordering in narrow energy bands is up to now extensively 
discussed. Despite the large number of publications on the topic, the magnetism of highly- 
correlated electronic systems described by the Hubbard model [1] remains at the center 
of attention [2—11]. Physically, in this case the picture of magnetism is characterized by 
the existence of local magnetic moments and differs substantially from the Stoner picture 
of a weak itinerant magnetism [12]. 

According to Nagaoka [13], in the limit of infinite Hubbard repulsion the ground state 
for simple lattices is a saturated ferromagnetic state for a low density 5 of current carriers 
(doubles or holes in an almost half- filled band). In particular, Nagaoka proved rigorously 
the existence of saturated ferromagnetic state for a single hole at U — *■ oo and found the 
instability of the spin-wave spectrum in the case with increasing S and decreasing U. 

Roth [14] applied a variational principle to this problem and obtained two critical 
concentrations. The first one, 8 C , corresponds to instability of saturated ferromagnetic 
state, and the second one, 8' c , to the transition from non-saturated ferromagnetic into 
paramagnetic state. For the simple cubic (sc) lattice, the values 5 C = 0.37 and 5' c = 



1 



0.64 were obtained. Next, the region of stability of the saturated ferromagnetic was 
investigated within various approximations in numerous works (see, e.g., Refs. [15-29]). 

An interpolational approach to description of magnetic ordering in narrow bands, 
which yields saturated ferromagnetism for small S and non-saturated one for large S, was 
developed in Refs. [30, 31] on the basis of dynamic magnetic susceptibility treatment. 
However, the critical concentrations themselves were not determined. 

Using high-temperature expansions in early papers [5,6,22,32,33] yielded non-stable 
results concerning stability of ferromagnetsm because of poor accuracy connected with 
slow convergence [18]. However, according to recent results [5,6], ferromagnetism also 
occurs near 5 = 0.3. 

It should be noted that the hole concentration 5 = 1/3 corresponds to the sign change 
of the chemical potential in the Hubbard-I approximation [1] in the case of symmetric 
conduction band, and instability of the paramagnetic state can occur at this point in 
some simple approximations [2]. 

Experimental data on Fei_ x Co :E S 2 [15], a system with strong correlations, give large 
values of 5 C (saturation ferromagnetism is preserved up to conduction electron concentra- 
tions n = 1 — 5 of order 0.2), but degeneracy effects in the conduction band appear to be 
important in this system. 

The approaches mentioned do not analyze as a rule in detail the structure of the 
one-particle excitation spectrum in the ferromagnetic phase of the Hubbard model. The 
simplest "Hubbard-I" approximation [1] for the electron spectrum corresponds to the 
zeroth order in the inverse nearest-neighbor number 1 / z ( "mean-field" approximation in 
the electron transfer). This approximation is quite non-satisfactory at describing ferro- 
magnetism (in particular, ferromagnetic solutions are absent, except for peculiar models 
of bare density of states). 

A consistent calculation of the one-particle Green's functions in the case of small 
5 (almost half-filled band) and low temperatures was performed in Refs. [34, 35] . The 
results demonstrated an important role of non-quasiparticle (incoherent) states in the 
density-of-states picture. Expressions for the one-particle Green's functions in a more 
wide region of 5 and T were obtained in Refs. [36]. 

A physically transparent mechanism of instability of saturated ferromagnetic state was 
treated in detail in the works by Edwards and Hertz [17]. This mechanism is connected 
with occurrence of spin-polaron states above the Fermi level of the current carriers. 

In this paper, the stability of the saturated ferromagnetic state as the current carrier 
concentration is raised is studied using the one-particle Green's functions of first order 
in 1/z and corresponding self-consistent approximations. This approach makes possible 
to construct a rather simple and physically transparent picture of the density of states 
in a saturated Hubbard ferromagnet. At the same time, the problem of description of 
non-saturated state is much more difficult, but our approach turns out to be successful 
too. 
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2 Calculation of the one-particle Green's functions 

We shall use the Hamiltonian for the Hubbard model in the limit of infinitely strong 
Coulomb repulsion in the many-electron X-operator representation [37], 

w = 5>*W (i) 

kcr 

where t k is the band energy, X^f is the Fourier transform of the Hubbard operators 
X" = \ia)(ij3\ (0 denotes holes and a = ±(f, |) denotes singly occupied states). 

It should be noted that in this problem of infinitely strong Coulomb interaction, 
a number of difficulties arise in connection with the non-Fermi excitation statistics. 
These difficulties occur both in the diagram technique [2] and in the equations-of-motion 
method [38]. In particular, it has been found [38] that in the expansion with respect 
to \jz the analytic properties of the retarded Green's functions were violated for the 
paramagnetic state. 

We shall calculate the retarded anticommutator Green's functions 

G ka (E) = ((XZ°\X°l)) E , %E>0. (2) 
We write down the equation of motion 

(E - t^)G ka (E) = (n + n a ) + r k(T (£), (3) 

where 

tka = t k (n + n a ), n a = (X° a ), a = 0, a, 

T ka (E) = ^2t k -q{{X^~ a X^ + + Xq <T )Xk° q |X°k)) J e, 

q 

with 5 A = A — (A). Decoupling the sequence of equations of motion at the first stage 
corresponds to the zeroth order in 1/z and is known as the Hubbard-I approximation. 
This may represented in the form 

Gl(E) = [F° a (E) - i k ]-\ F°(E) = - (4) 

When taking into account fluctuations we obtain 

r kCT (£) 



n + n a ) 

Following Ref. [36], we perform decoupling at the next stage to derive 

(E - * k _ q _ CT - auJ{{X^°X^\X%)) E = x7 a + ^- q -. 

+ - (*k- q - t k )n k _ q _ a )«X£° \X°l)) E , (5) 

(e - + x?)x?-Jix%))e = x-" 7 ^ + UX?\X%)) E ), (6) 

where 



are the correlation function for spin and charge densities 

n^ = (X%X?). 

The magnon frequencies in © are required to cut the logarithmic "Kondo" divergences 
which are connected with the Fermi functions. 

We write down the Green's function in the locator form 

GUE) = [F ka (E) - t k }-\ F ka (E) = (7) 
For the Green's function (JJJ) we have 

q 

+ J2t^X^- a G° k _ (lr7 (E), (8) 
q 

b ka (E) = £ + £4_ q n k _ q _ CT Gt q _ CT (£-^ q ). (9) 
q 

To simplify our equations, we introduce the magnon spectral function Kq(u) and average 
this in q, 

K^lu) = 6(u - cu q ) - K(u) = J2K n (uj). 

q 

This approximation is sufficient to obtain qualitatively valid results. Indeed, this ap- 
proximation (which is in spirit of the large-c? or large- z expansion) reproduces correctly 
the low-frequency behavior of spin fluctuations which is important near the Fermi level, 
following to Ref. [13] we have taken below D = 0.7S\t\. It should be noted that the 
choice of D influence weakly the critical concentration. We also neglect q-dependence of 
transverse and longutidinal spin correlation functions by taking the values averaged over 
the Brillouin zone to obtain 

Then a(E) and b(E) do not depend on k and can be expressed in terms of the bare 
electron density of states N (E), 

G k(T (£) = [F a (E) - t k )-\ F a (E) = Mfr, (10) 



a a (E) = n Q + n a + I K{u)^t^(n„ + n^ a )G°^_ a {E - au)du 

3 q 

+ £* q n_ CT (l-n_ CT )G qCT (£), (11) 
q 

K(E) = E + j K{u) ]T tln^ a G\_ a {E - aco)duj } (12) 
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G\_ a {E - auj) - G°_ a (E - acu,t) = [F a (E - au) - t}' 1 , F a (E) E 



n + 



a a (E) = n + n a + J K{uj) J N (t)t(n a + f_ (t))G°_„(E - ctlu, t)dudt 

+ J N (t)tn_ a (l - n_ a )G° a {E,t)dt, 
K{E) = E+ [k(oj) j ' N {t)t 2 f^{t)G°_ a (E -au,t)dudt, 

where f a (tq) = n qcr , N (t) = X)k5(£ ~ *k) is the bare density of states. The exact 
quasiparticle density of states is given by 

N a (E) = --$J2 G UE). (13) 
n k 

+00 

N a {E) = — 9 / G a (E,t)N {t)dt. 

7T J 



— 00 



Now we write down the self-consistent approximation with replacing in p ljl . (fT2|) the 
bare locator F®(E) by the exact locator F a (E), i.e. 

Gl(E) -> G ka (E). (14) 

In such an approach, large electron damping is present which smears the "Kondo" peak, 
so that including magnon frequencies does not change qualitatively the picture, but may 
be important for quantitatative results. We obtain 

GUE) = [F a {E)-^]-\ ^(£) = f^fj> (15) 

A a (E) = n + n a -\ ( K{u) Vt q (fi ff + n q _ CT )G q _ CT (£' - au)du 

n + n_ CT J q 

+ 7 X^n^l-n^GVCS), (16) 

^0 + n cr q 

B a {E) = E + f K{u)Y,tW-aG^ a {E - auj)duo. (17) 



A a (E) = n + n a + / K(u) f N {t)t{n a + /_ ff (t))G_ ff (E - ato, t)dtdu 

/ No(t)tn- a (l-n- a )G a (E,t)dt, 

n + n a J 



BJE) = E + / N (t)t 2 fJt)GJE - au, t)dtdu. 

Tin + U-rr J 



n + ti — (j 

It should be noted that another self-consistent approximation used in Refs. [38] leads 
to not quite satisfactory results because of violation of normalization condition for the 
density of states. As demonstrate our calculations, such a difficulty exists also for the 
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approximation (|10J) . but the violation is numerically small for 5 < 5 C (about 2%); note 
that introducing longitudinal fluctuations improves considerably the results in comparison 
with Ref. [39]. At the same time, the approximation (|THj) with the locator structure of 
the Green's function does not violate the analytical properties. 

The chemical potential /i is determined by the number of holes. By using the spectral 
representation for the Green's function (|2j) this is calculated as 

+oo 

6 = n = (X 00 ) = {X^Xf) = £<X^n> = — 9E / G ka (E)f(E)dE. (18) 

k 77 k 

It is important the Hubbard-I approximation is hardly satisfactory in the narrow-band 
ferromagnetism problem since it is difficult to formulate a reasonable criterion of magnetic 
ordering by direct using the expressions for one-electron Green's functions like Unlike 
the decoupling scheme by Hubbard [1], the many-electron X-operator approach clarifies 
the causes of this failure. In particular, one can see that the approximation (jlj) violates 
the kinematical requirements since it is impossible to satisfy at (S z ) ^ the relation 
(fT£j) for both spin projections a. Indeed, the quasiparticle pole for o =|, corresponding 
to a narrowed band and lying above the Fermi level of the holes, does not provide an 
adequate description of the energy spectrum and leads to the appearance of finite rij, 
i.e., the saturation ferromagnetism cannot be properly treated. However, the situation 
changes provided we use the expressions containing first-order l/^-corrections. Unlike 
the Hubbard-I approximation, the value of fi turns out to be weakly dependent of a for 
the approximation (fTU|) and independent of a for the the approximation (|15|). 
The magnetization is determined from the equation 

{S Z ) = \Y.^ (19) 



with 

n. 



+oo 

. = E^-W) = / N„(E)(1 - f{E))dE. 



In the leading approximation in 1/z for the one-particle occupation numbers, it is nec- 
cessary to use the Hubbard-I approximation, i.e., 

= {n + n a )f(t ka ), 

but the chemical potential should be already chosen from the Green's function ((7|). As 
opposed to Eq. the Green's functions ((Zj) contain terms with resolvents and have 
branch cuts which describe non-quasiparticle (incoherent) contributions to the density of 
states. It is the latter which ensure qualitative agreement with the sum rule (|T%|) for a =|. 
At the same time, there are no poles of the Green's function for this projection of the 
spin for small 5 above the Fermi level, i.e., the saturated ferromagnetic state is preserved. 

The full density of states can be also represented in terms of the exact resolvent 
R*(E) =£ k G kCT (£) as 

N a (E) = --%R a (E). 
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This quantity satisfies the equation 



R a {E) = R {F a (E)}, Ro(E) = £ 



1 



E-U 



N (t) 



1 



E-t 



dt. 



k - tk 

In the case of a saturated ferromagnetic state, the Green's function (JZj) takes the form 

-i 



G kT (£) 
G n (E) 



E ~E 

q 

x ( E 



1 + £ ^ 



E — t q n 

nk-q 



(20) 



, , v-^ (-B - £k)n k 
1 -" + ?^t 



t-k-q 



-1 



E ^k-q^k- 



(21) 



where n k = /(t k ). Note that, as opposed to the one-electron approach [17], the Green's 
function G^(E) does not reduce to the free electron Green's function even in the satu- 
rated ferromagnetic state, since fluctuations in the hole occupation number contribute to 
it. Neglecting the resolvent in Eq. (|2T)|) and the last term in the denominator of Eq. (|21jl. 
we obtain a somewhat different form of the Green's function in terms of the electron 
self-energy 



G kT (£) 



£-tk 

^l(E) = -(1 



1 



n 



E 



E — t k — T,ki(E) 



n k _ 



(22) 



E — t k _ q + UJ q 



This result corresponds to the Edwards-Hertz approximation in the limit U — > oo. Of 
course, these expressions work only in the saturated ferromagnetic state. The results 
(]22[) can be also obtained by using the expansion in the electron and magnon occupation 
numbers [40]. 



3 Results of calculations and discussion 

The 1/z corrections lead to a non-trivial structure of the total quasiparticle density of 
states. In the non-self-consistent approach the integral with the Fermi functions yields, 
similar to the Kondo problem, the logarithmic singularity 

E — — In — Ep\N(Ep). 

q - £k+q 

For very low S a significant logarithmic singularity exists only in the imaginary part of 
the Green's function, which corresponds to a finite jump in the density of states [34,35]. 
However, when 5 increases, it is necessary to take into account the resolvents in both the 
numerator and denominator of the Green's function, so that the real and imaginary parts 
are "mixed" and a logarithmic singularity appears in the density of states. When the 
magnon frequencies are included in the denominators of Eqs. (fTT]) and p2j1 . the singularity 
is spread out over the interval w max and the peak is smoothed out. In the self-consistent 
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approximations (fT6|) and fTTj) the form of Ni(E) approaches the bare density of states 
and the peak is completely smeared, even neglecting spin dynamics (Fig. Q), so that the 
latter plays no crucial role, although shifts somewhat the peak below the Fermi level. 

Near the critical concentration the peak in approximation (J7|) (but not in the Edwards- 
Hertz approximation) is again smeared (Fig. EJ), but this spreading out is no longer 
noticeable for 5 = 0.15. 

In the non-saturated state a spin-polaron pole occurs, so that quasiparticle states with 
o =1 occur above the Fermi level with 



The corresponding DOS picture is shown in Fig. |3J 

The critical concentrations 5 C for the loss of stability of saturated ferromagnetism, as 
calculated in different approximations considered, are listed in Table Q for a number of 
bare densities of states. In the case of fee lattices (where the bare density of states is 
asymmetric and has a logarithmic divergence on one edge) we have chosen the sign of 
the transfer integral for which the saturated ferromagnetism is stable at low 5 [13]. Note 
that it is necessary to use an equation for the chemical potential (|T8|) that is derived from 
the complete Green's functions (115)1 . (Using the Hubbard-I approximation here leads to 
a drop in S c of the order of 0.1.) 

It is clear from the Table ^ that the results are fairly stable and do not depend too 
strong on the form of the approximation. In particular, self-consistency changes them 
little, leading to a slight reduction in S c . The dependence on spin dynamics (magnon spec- 
trum), even in the non-selfconsistent approximation, is weaker (the critical concentration 
only varies in the third decimal place). At the same time spin dymamics is important 
for the description of the states near the Fermi level. It is interesting to note that re- 
sults of the Edwards-Hertz approximation (}2~2*|) are closer to those of the self-consistent 
approximation (|T3|l than of the non-self-consistent approximation ((7j). Unfortunately, in 
Ref. [17] only a crude estimate of S c was made by using the quadratic dispersion relation 
for the hole spectrum which yielded 5 C = 0.16. This approximation is not sufficient for 
quantitative calculations, as one can see from Tabled 

Unlike most other analytical approaches, our results for the one-particle Green's de- 
scribe formation of non-saturated ferromagnetism too. It is important that the account 
of longitudinal spin fluctuations (which were neglected in Ref. [39]) turns out to be 
important for obtaining the non-saturated solution and calculating the second critical 
concentration S' c 

The dependence of the saturation magnetization on the concentration of current car- 
riers for various bare DOS's is shown in Fig. 0] One can see that this dependence deviates 
from the linear one, (S z ) = (1 — n)/2, for 5 > 5 C . 

Let us perform a comparison of our results with other calculations. Generally, most 
calculations for a number of lattices yield the value of 5 C which is close to 0.3 (although 
the small value 5 C = 0.045 was obtained in Ref. [21] for sc lattice by a diagram approach, 
a close approach of Ref. [24] yields much larger values, e.g., 5 C = 0.25 for the quadratic 
lattice). At the same time, for the critical concentrations S' c the interval of values is 
broader and varies from 0.38 to 0.64. Our calculations yield the 5' c values which are 
considerably smaller than the results of the spin- wave approximation [14]. 
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Improved Gutzwiller method [26] yields for the sc lattice 5 C = 0.33, and using the 
t/U expansion [22] yields 5 C = 0.27. For the quadratic lattice the result of the variational 
approach [8] is 6 C = 0.251, and the result of Ref. [41] is 8' c = 0.38. The density matrix 
renormalization group approach [9] lead to the value 5 C = 0.22 and the rough estimate 
S' c ~ 0.40. Quantum Monte-Carlo (QMC) method for density matrix in Ref. [9] gives 
5 C = 0.22 and rough evaluation 5' c « 0.40. QMC in 2d case [42] gives 5' c « 0.40. Self- 
consistent spin density approximation (SDA) [43] lead to the results for simple cubic and 
bcc lattices 5 C < 0.32; the values of 8' c for the cubic lattices are given in Table [T] 

The values we have obtained can be compared with those in the limit of an infinite- 
dimensional space (it should be expected that our method of expanding in powers of 1/z 
is rather close to this approximation), for which 5 C = 0.42 (Ref. [26]) and S c = 0.33 
(Ref. [10]) have been obtained. At the same time, our approach is possible to reproduce 
the dependence of 5 C on the dimensionality of space and on the form of the bare density 
of states. 

Recently, 5 C has been obtained for a large number of lattices [3,28]. These results are 
also given in the Table ^ for comparison. It can be seen that in a number of cases our 
results agree better with a number of other calculations, especially for a square lattice. We 
note in this connection that a variational method has been used [4] to obtain a rigorous 
estimate of 5 C < 0.29 for a square lattice. Therefore, our results can be regarded as fairly 
reliable, even quantitatively. 

To conclude, we have obtained the density-of-states pictures in a Hubbard ferro- 
magnet with account of the "Kondo" scattering and spin-polaron contributions. Our 
approach yield a rather simple interpolational description of saturated and non-saturated 
ferromagnetism. One can expect that the results obtained will be useful for qualitative 
understanding of the ferromagnetism formation in narrow bands. 

The research described was supported in part by Grant No. 747.2003.2 (Support of 
Scientific Schools) from the Russian Basic Research Foundation. 

Table 1: Values of critical concentrations 5 C and 5' c for rectangular (ra) and semielliptic 
(se) bare density of state, square, simple cubic (sc), bcc, fee lattices. I, VII is non- 
self-consistent approximation (|7jl. II is Edwards-Hertz approximation (}2*2*|) . Ill, VIII is 
self-consistent approximations (with fluctuations) (fl5|) . IV is results of work [3], V, IX is 
results of work [43] VI is the result of Ref. [28] (variant of calculation RES0, for fee lattices 
instability not discovered). 



DOS 


S c 


K 


I 


II 


III 


IV 


V 


VI 


VII 


VIII 


IX 


ra 


0.276 


0.284 


0.301 








0.468 


0.488 




se 


0.258 


0.266 


0.290 








0.458 


0.477 




square 


0.244 


0.252 


0.275 


0.49 




0.4045 


0.449 


0.461 




sc 


0.233 


0.237 


0.261 


0.32 


< 0.32 


0.237 


0.427 


0.447 


0.66 


bcc 


0.217 


0.221 


0.247 


0.32 


< 0.32 


0.239 


0.414 


0.432 


0.48 


fee 


0.210 


0.217 


0.241 


0.62 






0.409 


0.427 


0.38 
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Figure 1: Density of states for the semielliptic DOS at concentration of carriers cur- 
rent S = 0.02. (a) line 1 (er =f) and line 2 (er =[) correspond to the non-self-consistent 
approximation (JJJ); line 1 (er =|) and line 3 (er =|) to the non-self-consistent approxima- 
tion with account of spin dynamics (fTU|) : (b) line 4 (cr =t) and line 5 (er =|) correspond 
to the self-consistent approximation f)15jl . and line 4 (cr =|) and line 6 (er =|) to the 
Edwards- Hertz approximation (|22|). 
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Figure 2: Density of states for the semielliptic DOS at concentration of carriers cur- 
rent S = 0.20. (a) line 1 (a =f) and line 2 (a =|) correspond to the non-self-consistent 
approximation with spin dynamics (jlO)) : (b) line 3 (a =|) and line 4 (a =!) correspond 
to the self-consistent approximation (|TK|) . and line 3 (cr =f) and line 5 (a =|) to the 
Edwards- Hertz approximation (|2*2*j) . 
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Figure 3: Density of states for the semielliptic DOS at concentration of carriers cur- 
rent S = 0.35. (a) line 1 (er =|) and line 2 (er =|) correspond to the non-self-consistent 
approximation with spin dynamics (b) line 3 (er =|) and line 4 (er =|) correspond 

to the self-consistent approximation Ijlfij) . 
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Figure 4: The dependence of the magnetization (S z ) on concentration of current carriers 5 
for a number of bare DOS's. Line 1 corresponds to rectangular DOS, line 2 to semielliptic 
DOS, lines 3, 4, 5, 6 to square, simple cubic (sc), bcc, and fee latices, respectively. 
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